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We investigated decoherence of a Josephson vortex quantum bit (qubit) in dissipative and noisy 
environment. As the Josephson vortex qubit ( JVQ) is fabricated by using a long Josephson junction 
(LJJ), we use the perturbed sine-Gordon equation to describe the phase dynamics representing a 
two-state system and estimate the effects of quasiparticle dissipation and weakly fluctuating critical 
and bias currents on the relaxation time Ti and on the dephasing time T^. We show that the critical 
current fluctuation does not contribute to dephasing of the qubit in the lowest order approxima- 
tion. Modeling the weak current variation from magnetic field fluctuations in the LJJ by using the 
Gaussian colored noise with long correlation time, we show that the coherence time Ti is limited 
by the low frequency current noise at very low temperatures. Also, we show that an ultra-long 
coherence time may be obtained from the JVQ by using experimentally accessible value of physical 
parameters. 

PACS numbers: 74.40.+k, 74.50.+r, 74.78.Na, 85.25.Cp 



I. INTRODUCTION 

Novel superconducting quantum bits (qubits), such as 
charge (i.e., Cooper-pair box)ji flux^ and phase qubits, 3 
are good candidates for quantum information process- 
ing because these can be manufactured, controlled, and 
scaled more easily. Tens of quantum oscillations had 
been observed in these qubits, but a low level of de- 
cay which yields thousands of coherence oscillations is 
essential for realization of quantum computation. The 
longest coherence time of 0.5 /is has been reported for the 
quantroniumi (i.e., hybrid charge-flux qubit), but longer 
time may still be necessary. As many decoherence sources 
reduce the quantum oscillations, the measured value of 
the coherence time for the superconducting qubits is sub- 
stantially shorter 5 than that predicted by the simplest 
models of decoherence and that needed for the operation 
of a quantum computer. This requirement of obtaining 
ultra-long coherence time in the presence of the interac- 
tion between the qubit system and noisy environment is 
one of major challenges. 

Understanding the mechanisms of decoherence became 
a focus of much attention recently, and it remains an im- 
portant challenge for the superconducting qubits. An 
ideal solution is to isolate the qubit system from un- 
controlled degree of freedom in its environment and in 
the device itself. However, it is difficult to isolate the 
qubit system completely from the decoherence sources. 
These sources include, but may not be limited to, back- 
ground charge fluctuations in the substrate, 6 fluctuations 
in the tunnel barrier which produce microscopic tun- 
neling resonance, 7 and fluctuating electromagnetic back- 
ground. Also, low frequency variation in the critical 
current is present in all superconducting qubits. One 
way to obtain the ultra-long coherence time is to use the 
Josephson vortex qubit (JVQ) since it msy be immune 



to these sources. 

Recently, JVQ has been proposed as a new supercon- 
ducting qubit. 8 This qubit has two important advantages 
over other superconducting qubits. First, coupling be- 
tween the qubit system and the decoherence sources is 
weak at very low temperatures. For example, a Joseph- 
son vortex (i.e., fluxon) in a uniform long Josephson junc- 
tion (LJJ) does not generate any radiation during its mo- 
tion and is almost decoupled from other electromagnetic 
excitations in the junction. Also the qubit is immune to 
fluctuations in the critical current. Consequently, quan- 
tum coherence can be maintained much longer than other 
qubits which are susceptible to these decoherence sources. 
Second, as the fluxon dynamics in LJJ is described by 
using the perturbed sine-Gordon equation^ the decoher- 
ence sources for the qubit may be easily identified. For 
example, two important sources are the quasiparticle dis- 
sipation and weak current noise. Hence, the coherence 
time may be estimated more easily, but the effect of these 
sources has not yet been estimated for the JVQ. 

The JVQ exploits the property of fluxon ji2*ii which be- 
haves as a quantum particle at ultra-low temperatures. 
The fluxon trapped in a controllable potential well in a 
single annular LJJ shows (i) energy quantization in the 
potential well and (ii) macroscopic quantum tunneling 
(MQT) from a metastable stated Also, it was shown 
that the two quantum states for the qubit can be created 
by using a heart-shaped^ annular LJJ and by trapping a 
fluxon in a magnetic field-controlled double-well poten- 
tial. These quantum states may also be created using 
a linear LJJ with two closely implanted defect sites in 
the insulator layer>i£ A Nb-AlO^-Nb junction may be 
used to fabricate the JVQ, as shown schematically in 
Fig. 1. The dimensions of the junction, compared to 
the Josephson length Xj, are L x 3> Xj and L y <C Xj. 
The separation distance £ between the defect sites must 




FIG. 1: A LJJ stack is shown schematically as alternating 
layers of superconductors (S) and insulator (7). L x and L y 
denote the dimensions in x— and y— direction, respectively. 
J B it) denotes the time dependent bias current density. The 
filled dashed circles represent the microresistors which behave 
as pinning centers for fluxon. 



be larger than the critical distance i 0l as discussed below. 
Also, the preparation of initial state and the read-out of 
the final qubit state may be performed by using classical 
circuits— 

As the effect of the decoherence sources in the LJJ is 
described by using the perturbed sine-Gordon equation, 
the coherence time T2, 
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for the JVQ may be estimated without making further 
assumptions about the nature of the qubit-environmcnt 
interaction. Here T\ is the relaxation time, and T^ is 
the dephasing time. Since T^ is most sensitive to the de- 
coherence sources, extending T^ for the superconducting 
qubits is important for quantum computing applications. 
We compute Ti by accounting for the two sources: (i) 
quasiparticle dissipation and (ii) weak current noise. It 
is noted that the JVQ may also couple to other sources, 
such as microwave and phonon radiation, but the effect 
of these sources is expected to be small as they lead to 
higher order— contribution in the perturbation expan- 
sion than that considered in the present paper. In the 
JVQ, the weak current noise represents the low frequency 
magnetic field and current fluctuations in the junction. 
We focus on slow fluctuations since the qubits suffer from 
the presence of strong low frequency noise sources. The 
effects of these two decoherence sources have been inves- 
tigated also for other superconducting qubits 1 — and were 
found to be important. 

In this paper, we show that the JVQ can yield ultra- 
long coherence time because it couples very weakly to 
noisy environment at low temperatures. Before proceed- 
ing further, we outline the main results, i) Starting from 
the perturbed sine-Gordon equation, we show that the 
critical current fluctuation does not couple to the JVQ 
within the lowest order approximation. Consequently, 
this fluctuation effect does not lead to decoherence of the 
qubit. ii) We show that Ti at ultra-low temperatures is 



determined by the low frequency current noise since the 
dissipation effect due to qubit- environment coupling is 
exponentially small, iii) Using experimentally obtained 
physical parameters, we show that the effect of this cur- 
rent noise on decoherence is weak in the JVQ. This weak 
coupling between the JVQ and current noise leads to the 
coherence time of several tens of microseconds for the 
JVQ. 

We outline the remainder of the paper. In Sec. II, we 
express the phase dynamics of LJJ in the collective coor- 
dinate representation and transform the perturbed sine- 
Gordon equation onto the double- well potential problem. 
In Sec. Ill, we obtain the two-state system, described by 
the spin-boson model with low frequency current noise. 
Here, the quasiparticle dissipation is described by us- 
ing Ohmic environment, and the effects of low frequency 
noise in LJJ is described by using the fluctuating weak 
bias current. In Sec. IV, we derive the relaxation time 
(Ti) and the dephasing time (T^) due to these two deco- 
herence sources. In Sec. V, we estimate numerically the 
coherence time (T2) by using the experimentally accessi- 
ble parameters for LJJ and compare it to that obtained 
for other superconducting qubits. Finally, in Sec. VI, we 
summarize the result and conclude. 



II. LONG JOSEPHSON JUNCTION QUBIT 

In this section, we discuss how the LJJ may be used 
to obtain a JVQ by starting from the perturbed sine- 
Gordon equation. First, a double- well potential for the 
fluxon needs to be created in the LJJ to obtain the two 
quantum states of the JVQ. Several approaches are used 
to accomplish this. Each of these approaches yields a 
slightly different form of the potential function. In this 
paper, we consider the approach of implanting two closely 
spaced microresistors in the insulator (I) layer. When the 
fluxon does not have enough kinetic energy, the microre- 
sistor attracts the fluxon and traps it at the defect site. 
The effects of quasiparticle dissipation and low frequency 
weak current noise may be examined by starting with 



d 2 <p d 2 p . 
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where x and t are the dimcnsionless coordinates in 
units of A j and u>~ , respectively. Here io p denotes the 
plasma frequency The dynamic variable ip represents 
the difference between the phase of order parameter for 
the superconductor (S) layers. The perturbation term 
T = P(dtp/dt) - (3 s (d 3 <p/dtd 2 x) + f(t) + 8J c (t)smcp - 
'^2 i CiS(x — x°) sin p includes the effects due to quasipar- 
ticle dissipation (/? and (3 S ), bias current (/ = J B /J c ), 
critical current fluctuation (SJ c (t) — SJ c (t)/J c ), and mi- 
croresistors (ei = (J c — J' c )lb/J c Xj). We note that the 
critical current J c (t) may be expressed as the sum of 
uniform (J c ) and weak fluctuation parts (5J c (t)): J c (t) = 
J c + SJ c (t). Here x°, J B , J' c , and lb (<§^ A,/) denote, re- 
spectively, the position of inhomogeneity in the / layer, 



the bias current density, the modified critical current den- 
sity at the defect site, and the length of the LJ J in which 
J c is modified. In the discussion below, we neglect the 
term (3 s (d 3 cp/dtd 2 x) due to the quasiparticle (surface) 
current along the junction layer, for simplicity, since both 
(3{dip/dt) and f3 s {d 3 ip/dtd 2 x) terms yield similar dissipa- 
tion effects. 

The perturbation term T in Eq. @ is small, and con- 
sequently, it does not change the form of the kink within 
the framework of the lowest approximation. 9 We describe 
the motion of the fluxon in terms of the center coordi- 
nates q(t), which are obtained by neglecting T. In the 
absence of the perturbation terms [T = 0), the fluxon so- 
lution to Eq. J5J in the non-relativistic limit (i.e., !)«1) 
is given by 



tp(x, t) fa 4 tan 



,-y(v)[x-q(t)] 



(3) 



where r y~ 1 (v) = yl — v 2 , q(t) — vt denotes the center 
coordinate of the fluxon, and v is the fluxon speed in 
units of Swihart velocity c. q(t) is also known as collec- 
tive coordinate and represents a dynamical variable. We 
note that this solution represents propagation of nonlin- 
ear wave as a ballistic particle and that the perturbation 
terms in J- only affect dynamics of the center coordinates. 
Applying the kink solution of Eq. © to Eq. (0) for 
carrying out the classical perturbation theory within the 
framework of lowest order approximation, 9 we obtain the 
equation of motion for the center coordinate q(t) in the 
nonrclativistic limit as 



M 



d 2 q(t) 
dt 2 



(3M 



dq(t) , dV(q) 
dq 



dt 



= 



(4) 



Here M — 8 denotes the rest mass of the fluxon which 
is obtained by inserting the waveform of Eq. J3J into 
the Hamiltonian corresponding to the unperturbed sine- 
Gordon equation (i.e., Eq. J2J with T = 0). 9 V(q) is 
the effective potential for the fluxon due to the non- 
dissipative perturbation terms in T . 

The phase dynamics in the center coordinate may be 
seen easily from the Euclidean Lagrangian, £ = £ + Cp 
where C and Cp describes the unperturbed phase dy- 
namic of LJJ and the perturbation contribution, re- 
spectively. The unperturbed part of the Lagrangian 
C is given by C = J(dx/2)[(dtp/dT) 2 + (d<p/dx) 2 + 
2(1 — cos ip)]. The perturbation part of the Lagrangian 
Cp = C n d + Cd can be expressed as the sum of the non- 
dissipative part (C n d) due to the bias current, critical 
current fluctuation and inhomogeneities, and the quasi- 
particle dissipation part (Cd)- The non-dissipative con- 
tribution may be expressed as C n d = Cbias + Cgj c + C P i n . 
Here Cbias = J dxfip, C j c = SJ c (t) J dx(l — cosip), and 



2 

Cpi n — £ / dx ei6(x - x°)(l - cos ip) 

i=l 



(5) 



are the Lagrangian for the bias current, critical current 
fluctuation, and two defect sites, respectively. Following 



Caldeira and Leggett^ we account for the quasiparticle 
dissipation (i.e., (3) by representing the environment as a 
heat bath. The heat bath may be described as a reservoir 
of harmonic oscillators with generalized momenta Pi and 
coordinates Qi. The dissipation effect due to coupling 
between the phase variables ip and the heat bath is de- 
scribed as 



C d 



/<**£ 



2m,- 



Qi 



Cj<P 
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Here, the heat bath parameters m,-, u>i and q characterize 
the reservoir's spectral function Jp(oj), which is written 
as 



I t— 1 rriiUJi 



6(lu — oji) = (5u> 



(7) 



This spectral function reproduces the dissipation term 
in Eq. Q when the heat bath degrees of freedom are 
integrated out. 

We describe the fluxon dynamics by using semiclas- 
sical theory as usually doneii and by reexpressing the 
partition function, Z = J T>[ip] exp{— S* [</?]}, with S[tp] = 
J drC, in terms of the collective coordinates q(t) as 
Z = J T>[q]cxp{—S[q]}. We take the perturbation ex- 
pansion in terms of /3, /, 5J C and e,, assuming that all of 
these parameters are small. The lowest order contribu- 
tion from this expansion is obtained by substituting the 
soliton solution of Eq. © to the action S[ip] since the 
perturbation term T does not modify the soliton wave- 
form in the lowest orderiSiil In this center coordinate 
representation, the bias current contribution to the ac- 
tion (i.e., SuasYp] = f dr Cbias) yields 10 



Sbias[q] = S bias [tp(x - q)} - S b ias[<p{x)} =- dT(2wfq) . 

Here the constant Sbiasipfa)] is subtracted since we have 
chosen the origin of potential energy for the center coor- 
dinate at q — 0. On the other hand, the critical current 
fluctuation contribution (i.e., Ssj c — J drCsj a ) becomes 



SsjM = Ssj e [<p(x - q)} = Ssj c [ip(x)} 



(9) 



indicating that S$j c [q] is independent of q. This sug- 
gests that the critical current fluctuations do not modify 
the fluxon potential. Within the lowest order approxima- 
tion, the low frequency noise corresponding to the criti- 
cal current fluctuation does not couple to the JVQ, and it 
does not contribute to decoherence. Also, other perturba- 
tion contributions Spin = J drC V i n and Sd = / drCd can 
be expressed in the q representation. Combining these 
perturbation contributions, we may express the partition 
function as Z = j T)[q(T)]exp{—S e ff[q(T)]}, where the 
effective action S e ff [q] is given by 



S eff[q] 



dr 



-Mq 2 + V(q) 



+ ^-JdrJdr%(r - r')[q(r) - q{r')f 



(10) 
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FIG. 2: A double-well potential due to two microresistors is 
schematically shown for (a) ei = 62 = e (symmetric) and (b) 
£i 7^ £2 (asymmetric). The potential barrier and the oscilla- 
tion frequency at the stable minimum are denoted by V and 
uj , respectively. Here f(t) is set to zero. 



The quasiparticle dissipation effect (i.e., /?) at the finite 
temperature T is described by the kernel K.(t) 

1C(t) = -J o ** W Sinh(h;/2T) ■ (ii) 

Here we set 7i = fes = c = 1 for convenience. The poten- 
tial function V(q) for the fiuxon in the collective coordi- 
nates is given by 



V(q) = -2nf(t)q- 



2 £l 



2e 2 



cosh [q — |) cosh [q + |) 



(12) 



where I is the separation distance between two defect 
sites, as shown in Fig. 1. The fiuxon potential V(q) 
includes the potential tilting effect of the bias current 
(/) and the pinning effect (ei) of the two defect sites. 

The potential function V(q) of Eq. (JEJ for /(*) = 
has two noteworthy features: (i) finite number of bound 
states, and (ii) double-well structure. For physical values 
of e, , at most, several states may be trapped by the fiuxon 
potential. This can be seen easily from the energy eigen- 
state of the trapped fiuxon via a single microresistor, 18 
which is given by 

£ n = ~ [-(1 + 2n) + \/l + 128e f (13) 



where n = 0, 1, 2, • • •. For e = 0.27, only n —0, 1, and 
2 states, corresponding to the eigenstate energy £ =- 
0.383, £\ =-0.136, and £2 =-0.0142, respectively, are 
bounded by the potential. Also, the double-well struc- 
ture can be seen easily by setting e\ = e 2 = e (i.e., sym- 
metric double- well) (see Fig. 2(a)) and by expanding the 
function V(q) about the critical separation distance £ . 
We note that a small asymmetry (or bias) of e ~ 8q Se 
may be easily introduced, as shown in Fig. 2(b), since the 
critical current J' c at the each microresistor is slightly dif- 
ferent. This yields a small variation in &j between the two 
defect sites (i.e., |ei — e 2 | = Se and Se <C 1). The symmet- 
ric potential V(q) shown in Fig. 2(a) has the single-well 
structure for £ < £ = ln[(\/3+ 1)/(V3- 1)] ~ 1-317, but 
it has the double- well structure for £ > £ , For £ = £ + a, 
with s«l, the fiuxon potential V(q) may be expanded 
around q = to obtain 



V(q) - V(0) 



16e 2 32e 4 
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aq 



27 



(14) 



The potential function of Eq. (|14f> has the stable states 
at q = ±q /2 with q — 3[a/v / 3] 1 ^ 2 - The barrier height 
V between two stable states is V — 2ea 2 and the fre- 
quency lj of small oscillation around the stable mini- 
mum is u = y/(d 2 V(q)/dq 2 )/M = {8ea/SV3) 1/2 , indi- 
cating that lo ^> V . This suggests that the JVQ can- 
not be obtained by placing two microresistors too closely 
(i.e., £ rs £ ) since the potential barrier V may not be 
strong enough to localize the fiuxon to either well. Hence 
a larger separation distance £ is needed to localize the 
bound states in either potential well. 

The barrier potential V must be larger than to in or- 
der to obtain a localized ground state in either well with- 
out mixing it with the excited state of the system. As £ 
increases from £ , the value of both V D and lo increases, 
but this increase depends on e . As £ — > 00, V approaches 
2e while cu approaches y/e/2. This indicates that, when 
e is less than the critical value e c (i.e., e < e c ), V re- 
mains smaller than io for all £. However, when e > e c , 
V becomes larger than lu as £ is increased, as shown 
schematically in Fig. 3. We estimate e c = 0.125, assum- 
ing that uj (e c ) — V (e c ) at £ = 00. We will consider 
e > e c in the discussion below, so that V > lu . 

For the fiuxon localized at either left or right side of 
the symmetric double- well shown in Fig. 2(a), the energy 
of the ground state is degenerate. We use eigenstates \R) 
and \L) of the operator a z with eigenvalues +1 and -1 
to represent the right-localized and left-localized state, 
respectively. These two states are exploited in the JVQ. 
To do this, we need to ensure that the tunneling rate 
between the two wells does not mix the ground state 
with the excited states. 

The fiuxon in the ground state of the symmetric 
double-well potential can tunnel from the left side to the 
right side (and vice versa) . This MQT yields splitting of 
two degenerate fiuxon ground states. Within the semi- 
classical WKB approximation^ the tunneling rate A is 
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FIG. 3: Dependance of the barrier potential V (dashed line) 
and the frequency of small oscillations u> (solid line) for the 
symmetric double-well potential on the separation distance £ 
between the two defect sites is schematically illustrated for a) 
e < e c and b) e > e c . The critical pinning strength e c is 0.125. 



given by 



A =.4(0) 



-8(0) 



where 



(&4foY ?l7 d * 



.4(0) 



B(0) = 4 / dq JV(q) , 



(15) 

, (16) 
(17) 



and V(q) = V(q) — V(±q /2) denotes the potential en- 
ergy measured from the bottom of the well. The MQT 
in the real space represents particle-like collective exci- 
tation, reflecting the behavior of the fluxon as a quan- 
tum particle. The computed tunneling rate, using the 
WKB approximation, yields good agreement with the 
quantum result when many states are bounded by the 
double-well potential, but this agreement is poor when 
only the ground state is bounded. In Fig. 4, we com- 
pare the result of the semiclassical WKB calculation us- 
ing Eqs. (|15|) - (|17fl and the quantum mechanical calcula- 
tion to illustrate this difference. The difference between 
these two results is noteworthy: the WKB calculation 
overestimates A. This difference is large when £ is small 
(i.e., small V ) but decreases with increasing I (i.e., in- 
creasing V ). We will use the quantum result for A in the 
discussion below since the coherence time depends on A. 
Numerical solution of the bound state energy for the 
potential of Eq. I|12fl (with /(£) = 0) indicates that the 
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FIG. 4: The splitting A of the ground state for e = 0.27 is 
plotted to compare the results obtained from the WKB ap- 
proximation and that obtained from the quantum mechanical 
calculation. 




FIG. 5: The diagram in the (£, e) parameter space illustrating 
that bound fluxon states localized in either side of the double- 
well potential. The shaded areas represent the region with 
large tunneling rate. 



ground state is localized at the either side of the double- 
well for only limited value of £, as shown in Fig. 5. The 
lower and upper shaded areas in Fig. 5 represent the 
regions in the (e, £) parameter space where the tunneling 
rate between the two degenerate ground states and the 
excited states, respectively, are large so that these states 
cannot be localized in cither well. For a fixed e, the 
number of localized states in either well increases with 
£. This indicates that the separation distance £ and the 
pinning strength e may be chosen so that only the ground 
state is localized in either well. The 'parameters which 
yield localization of only the ground state (i.e., between 
two shaded regions) may be ideal for obtaining the JVQ. 



III. SPIN-BOSON MODEL 

In this section, we describe the interaction between the 
JVQ and noisy environment. We proceed by describing 
the fluxon dynamics of Eq. (|10|) in terms of the well- 
known spin-boson model. This may be carried out by 
using the two-dimensional Hilbert space spanned by the 
two degenerate ground states: the fluxon localized at the 
left well (i.e., \L) ) and at the right well (i.e., \R) ). Fol- 
lowing earlier studies, we consider the parameter regime 
of V > 0J o ^S> A, e, T and include the effects of quasipar- 
ticle dissipation and fluctuating weak bias current in the 
spin-boson modeli22i2i The Hamiltonian for this model is 
written as 

H = H S + H SB + H B . (18) 

The spin (S) Hamiltonian Hs, 

ffs = -~Aa„-i[e + /(«)]*,, (19) 

describes the two-state qubit system, which is obtained 
from the double- well potential of Eq. i|12[l . Here A is the 
tunneling rate between the two wells. The Pauli opera- 
tors, <t z and a x , in Eq. Q19[) represent 



and 



& Z = \R)(R\-\L)(L\ 



a x = \R)(L\ + \L)(R\ , 



(20) 



(21) 



respectively. The Hamiltonian Hs also accounts for a 
modification of the simple two-state system by a small 
asymmetry in the potential due to slight variation in the 
pinning strength of the microresistors (i.e., ei ^ £2) and 
by fluctuating bias current (i.e., f(t)). The bias current 
density f(t) = q f(t), representing the driving force for 
the fluxon, consists of two parts 



/(*) = fo + Sf(t) 



(22) 



where f and Sf(t) denote the homogeneous and ran- 
domly fluctuating weak bias current components, respec- 
tively. Here 6f(t) accounts for the current noise in the 
JVQ. The Ohmic environment^ which accounts for the 
quasiparticle dissipation, is described by the bath (B) 
Hamiltonian Hb, 



H B = 



1 N 

IE 



2 f—f \ in, 



p 2 



rmwlQl 



(23) 



The interaction between the qubit system and the dis- 
sipative environment is described by the spin-bath (SB) 
Hamiltonian Hsb, 



Hsb = -Vz-^y^CjQ 



N 



(24) 



It is noted that the spin-boson model of Eqs. I|18fl . l|19fl . 
(|23[) . and (|24|l neglects the contributions from the excited 
states that are bounded by the potential well. Hence, 
thermally activated leakage j22 which may also contribute 
to decoherence at finite T, is not accounted in this work. 
However, we may, safely, assume that this contribution is 
negligible at ultra-low temperatures. Consequently, the 
weakly fluctuating bias current (i.e., current noise) at low 
frequency is the dominant source for dephasing at these 
temperatures. 

We now discuss the time dependent bias current fit) 
of Eq. Ij22(l . We set the externally applied homoge- 
neous component of the bias current to zero (i.e., f = 0) 
since it yields unwanted asymmetry in the double-well 
potential (see Fig. 2(b)) for the fluxon. The weak bias 
current fluctuation (Sf(t)), representing random force in 
the LJJ due to nonequilibrium states, yields small time- 
dependent asymmetry in the double-well potential. This 
asymmetry may be made small but cannot be turned off 
completely as it arises from the noise-producing environ- 
ment. This current noise leads the basis states {\L), \R)} 
to fluctuate weakly, as shown schematically in Fig. 6. 
However, this effect on the basis state in the JVQ is ex- 
pected to be smaller than that in other superconducting 
qubits since the bias current is not used to control the 
qubit state. It is noted that the effect of noise in the bias 
current for f ^ has been investigated for the phase 
qubrbSii and for the charge qubit iSi In these qubits, the 
bias current (i.e., f ) is used to control the qubit. Noise 
in the bias current (i.e., 8f) affects the coherence time 
of the qubit since it leads to the fluctuation of the qubit 
state. The bias current noise leads to phase noise for the 
phase qubi1>2£ and radiation noise for the charge quibit^ 

To model the effect of this random force more realis- 
tically in the JVQ, we describe the bias current fluctua- 
tion as Gaussian colored noise with non-zero characteris- 
tic correlation time r„. The current noise 8 fit) has two 
main effects on the dynamics of qubit density matrix: 
i) it leads to transition between two energy eigenstates, 
and ii) it suppresses coherence between the eigenstates 
by contributing to pure dephasing. The Gaussian col- 
ored noise can be used to account for the noise spec- 
trum with pronounced frequency dependence, such as 
Lorentzian noise and low frequency asymmetrical mag- 
netic field fluctuations^ in the tunnel junction. The 
characteristics of the bias current are described as 



and 



<<*/(*)> = 



(8 fit) Sf(t')) = n\ e -l*-*'l/rn 



(25) 



(26) 



Here (• • •) denotes average over different realizations of 
the fluctuating current, and n is the typical noise am- 
plitude. We note that, as r„ — ► 0, the colored noise of 
Eq. H2(j|) becomes the white noise, which is characterized 
by the correlation function {8 fit) 8f{t')) = n^Sif. — t'). 
The correlation function of Eq. 126f> indicates that the 




FIG. 6: A schematic diagram illustrating the effect of colored 
noise on the bounded ground state of Josephson vortex in the 
symmetric double-well potential. 



effects of the current noise for t <C r n differ from those 
for t ^> r n . For Kr„, decay of coherence arises from 
averaging over the distribution of current noise since the 
fluctuations appear static. For t 3> T n , on the other 
hand, decay of coherence is expected to be exponential 
since the fluctuating bias current behaves as white noise. 
The crossover behavior occurs at t ~ r n . The spectral 
density of the bias current noise can be taken as 



o ,, 



M = 



1 + {uT n y 



(27) 



The Lorentzian spectrum of Eq. (|27() , characterizing tele- 
graph (diachotomous) noise, was observed in intrinsic 
LJJ. 26 We note that the noise spectrum in a small tunnel 
junction is described by the Lorentzian function of Eq. 
(|27|l . but the 1/w-like noise spectrum in a larger junc- 
tion may be obtained as a result of several superimposed 
Lorentzian features^ In the discussion below, we make 
few assumptions about the correlation function of Eq. 
(1201 : the fluctuation is weak and has small characteristic 
anrplitude (i.e., n <C A), but it has a long correlation 
time (i.e., r n 3> 1/A). Also, we assume that the temper- 
ature (T) of the bias current producing environment is 
larger than the cutoff frequency of l/r n (i.e., T ^S> l/r n ). 



IV. DECOHERENCE DUE TO FLUCTUATING 

WEAK BIAS CURRENT IN OHMIC 

ENVIRONMENT 

We now discuss the effect of dissipation and noisy envi- 
ronment on the coherence time of the JVQ by using well- 
established formalism. Here the coherence time repre- 
sents the time scale for decay of macroscopic quantum co- 
herence (MQC) between the ground states in the double- 
well potential. Here MQC is due to quantum tunneling 
of the fluxon which leads to coherent oscillations. This 
MQC is suppressed by the two decoherence sources since 
the interaction between the qubit system and its envi- 
ronment can easily destroy the phase coherence between 
two states. In estimating T2, we follow the standard the- 
oretical approach of using the Bloch-Redfreld theory and 
making lowest order Born approximation. The effects of 



these two sources may be characterized as follows. The 
Ohmic environment yields the finite relaxation time (Xi) 
and dephasing time (T?). However, the fluctuating weak 

bias current modifies T? without significantly changing 
T\. We estimate the effects of the bias current noise, 
restricting our consideration to t > t„ since the bias 
current fluctuation appears as 5- function correlated (i.e., 
white noise), and the quantum coherence decays expo- 
nentially. Hence Ta, may be expressed simply as 



1 



1 



1 



rpB rpnoise 

J 1 4> 



(28) 



where TV; mse is the dephasing time due to weak bias 
current noise. This indicates that the divergence in 
T¥ at ultra-low temperature may be cut off by TV; mse . 
The contributions from the higher-order Born correc- 
tion and non-Mar kovian effect, which are not included in 
the present work, are small but yield power-law decayS 
These contributions may also cut off the diverging T? 
due to significantly reduced interaction between the qubit 
system and environment at ultra-low temperatures. 

The decay of coherent oscillations is estimated by con- 
sidering the generalized master equation for the system's 
density matrix ps(t), 



dpsjt) 
dt 



= -i[H s ,p s (t)}-i f dt'X s (t-t')p s (t') , (29) 
Jo 



and by assuming that the time dependence of the de- 
coherence source is weak. Here ps(t) = Trspit), P — 
Ps® Pb-i the kernel Sg(t) is the self-energy operator 



£ s (i) = -m B H SB e- QHt H S B Pb 



(30) 



Hsb is the Liouvillian superoperator defined by TLsb P = 
[HsbiP], and Q = 1 — pb^tb is the projection superop- 
erator. Here pb is the bath density matrix. Since the 
studies^ indicate that both Bloch-Redfield theory and 
path integral theory yield equivalent results, we employ 
the former approach for convenience. 

The generalized master equation within the Born ap- 
proximation is obtained by using the fact that the cou- 
pling between the qubit system and environment, as de- 
scribed by Jp{ui)/uj ex (3, is small at low temperatures. 
We follow the Redfield theory^ and make a system- 
atic perturbation expansion of the kernel S5 in pow- 
ers of the system-bath coupling ((3). We retain only 
the lowest order terms in this expansion. Replacing 
e -iQHt _^ e —i{H s +H B )t anc i keeping the expansion of the 

kernel £5 up to the second order in Hsb, we obtain 



.(2) 



T,%>(t) = -iTr B H S Be 



-i{H s +n B )t 



T~(-sb Pb 



(31) 



Further simplification of Eq. I|29|) may be made by as- 
suming Markov system dynamics 



Ps(t-T ) ~ e 



iHst 



PS it) 



(32) 



in which the temporal correlation time r in the dissi- 
pative environment is very short due to very short-lived 



system-bath interactions, and the bath correlation func- 
tion decays to zero at a very short-time. 

Bloch-Redfield equation: In examining the decoherence 
effects, the Hilbert space spanned by the ground states 
of the two wells (Fig. 2(a)) is not convenient since the 
spin Hamiltonian of Eq. I|19(l is not diagonal in the basis 
{|£},|-R)}. We represent the two-state system in new 
basis {|0), |1)} given by 



|0) = -\L)sin6 + \R)cos9 
|1) = |L)cos6»+ | i?) sin 6» . 



(33) 
(34) 



where sinA = ^/(fi - e)/ti/V2, cos 9 = y/(Q + e)/£l/y/2, 
and Q — \Jl 2 + A 2 . In this new basis {|0), |1)}, we esti- 
mate T2 by making the Born-Markov approximation and 
by obtaining the Bloch-Redfield equations. Taking ma- 
trix elements in the eigenbasis \n) of Hs (i.e., |0) and 
|1)), we may write the Redficld equations as 

dpnmit) 



dt 



= -iE„ 



P„ 



nm Hum 



,(t) -J^Knmkl p S kl(t) , (35) 
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where p„ 



s _ 



\n\ps\m) 



K,, 



1 ki 



is the Redfield tensor, 



E n m = E n — E m , and E n is the eigenstate energy of 
Hs (i.e., Hs\n) = E n \n)). In the absence of the fluc- 
tuating bias current (Sf), the eigenstate energies in the 
diagonal basis, representing the ground state splitting, 
are expressed as 



E = -Iv/e 2 + A2, 



E x = +i v / e 2 + A2 



(36) 
(37) 



In the presence of the low frequency bias current fluctu- 
ations, these energies fluctuate slowly as shown schemat- 
ically in Fig. 6. We assume that the eigenstate energies 
are almost constant in the time scale relevant for the 
evolution of the density matrix. 

The Redficld tensor lZ nm ki is defined by 



K 



nmkl 



dm B (n\[H SB (t), [H SB (0), \k(t))(l(t)\p B ]]\m) (38) 



where the spin-bath Hamiltonian Hsb and qubit system 
eigenstate |fc) in the interaction picture are written as 



H SB {t) = e ^s+H B )t HsBe - t (Hs+H B )t ( (39) 

|jfe(t)) = e UHs \k)=e ttEk \k) , (40) 

respectively. The Redfield tensor may be expressed as 

1? —a Y^r(+) -la Y^r(~) _ p(+) _ p(-) 

K-nrnkl — Olm ^ 1 n rrk "+" ° nk 2-^t Irrm l Irnnk l Irnnk 
r r 

( 41 ) 

by evaluating the commutators in Eq. I)38|l. Here 



r 



(+) 

Irnnk 



dt e-^Tr B HtJ?(t)H^(0) p B ,(42) 



Imnk 



f 

Jo 



dt e- u ^Tr B H^(Q)H^(t) p B ,(43) 



and RS£(t) 



a itH t 



H SB e~ 



it H B 



The relation 



(T (+) 1 

V ImnkJ 



r 



(-) 



(44) 



may be used to write the Redfield tensor in terms of only 
the complex F im ^ fe tensor 



-.(+) 



2q 2 (l\a z \m)(n\a z \k) / — e 



dt 



-iE„ k t 



dujJfjiuj) 



2tt 



coth — — cos tot 
2T 



i sin tot 



(45) 



with the spectral function Jp(ui) of Eq. Q representing 
Ohmic environment. 

The dynamics of the two-state system may be de- 
scribed by using a 2-by-2 density matrix which is written 
in the Bloch vector form (i.e., three real variables). The 
Bloch vector p is written as 



P01 



P = Tr((?ps) = I i(p? 




(46) 



where a = (cr x , a y , cr z ) represents the vector composed of 
the three Pauli matrices. The Bloch vector of Eq. I|46|) 
may be combined with the Redfield equation of Eq. I|35|) 
to obtain the Bloch-Redfield equation 



dp 
-dt =eXP 



Rp + Po 



(47) 



where e = (0,0, Eqi) t , the relaxation matrix R is given 
by 



n[ 



R 



and 



^011 



0101 



0110 



2K'r 



0001 



Po 



"^•0101 — 'MH10 


^•0100 — 7^0111 


^-0101 — 'MlllO 


^•0111 — ^-oioo 


2^-ooni 


^•0000 + "-1111 




(48) 


(^•01 11 + ^-0100) \ 




^•0100 + ^0111 


(49) 


(^-0000 ~ ^1111) / 





Here Kmkl and K 



nmkl are the rea ^ and the imaginary 
part of the Redfield tensor, respectively. 

The Bloch-Redfield equation of Eq. I|47(l may be sim- 
plified within the secular approximation. Within this 
random phase type approximation which corresponds to 
retaining only the terms TZ nm ki with the indices n — m = 
k — I, the Redfield tensor simplifies to R ~ R se c and 



R, 



-K. 




n' 



0101 



0101 




(50) 



n' r 



0000 



n[ 



in , 



This approximation is valid for the spin-boson model of 
Eq. (JTBJ). We note that E 01 = ft while Tl nm ki < OQ3), 
when n — m^k — l. Since f2 3> P at ultra-low tempera- 
tures, Eqi > TZnmki- 



Ohmic environment: We solve the Bloch-Redfield 
equation of Eq l|47|l within the secular approximation to 
obtain the relaxation and dephasing time due to the en- 
vironment. The decay of the diagonal element of the 
qubit's reduced density matrix is written as 

p s z (t) = /of (0) e-*/ Tl . (51) 

This yields the relaxation time (T±) which is given by 



rr K 



oooo 



^i m -2Re(r + ) +rS ( ] 1 ). (52) 



Evaluating the tensors r ollo and r^goi by using Eq. Q45JI . 
we obtain the relaxation time T± as 



1 4A 2 



n 



q n B coth — 
Ti n HoH 2T 



(53) 



This is consistent with the result from the bounce 
solution^ The decay of the off-diagonal element of the 
reduced density matrix (either Pq X or pf ) may be ex- 
pressed as 

pS^t) = p*(0)I(t) e- f/T " e +^+<ioi]t . (54) 

This reduced density matrix includes the effects from 
both the Ohmic environment (i.e., XrP) and the fluctu- 
ating weak bias current (i.e., 1(t)). The coherence time 
T 2 B due to Ohmic environment is obtained as 



rpB 

1 2 



^A/^ L QrrO T 2-^1 lrrl 1100 L 1 



100 ) 



2Ti T^ 



(55) 



From Eq. (|55() . it is straightforward to obtain the de- 
phasing time (T?) due to the Ohmic environment as 



16e 2 

n 2 



q 2 oPT. 



(56) 



This indicates that T? diverges as either e or j3 vanishes. 
We note that e is a temperature independent parameter, 
but (3 becomes exponentially small at ultra-low temper- 
atures, yielding strong divergence in T^ . 

Fluctuating weak bias current: The dephasing time 
rpnoise < j me ^ Q ^- ne wea ]j bj as current noise may cut off 

the divergent T? at ultra-low temperatures. We esti- 
mate T£ mse from I(t) of Eq. (52J. The suppression 
factor X(t) in the off-diagonal clement of reduced density 
matrix p^ x represents the decay 3 ^ of coherence due to the 
bias current noise. This suppression factor 



I(t) =cxp\ ±i / dt' 



ef(t') , A 2 f 2 (t>) 



n 



2Q 3 



(57) 



accounts for the accumulation of the noise induced phase 
between two instantaneous energy eigenstates |0) and |1), 



due to long correlation time t„. We estimate T? mse by 
averaging T(t) over the realization of the fluctuating bias 
current and obtain 



I{t) = ( exp\±i alt' 



sf(t>) , A 2 f 2 (t<) 



O 



2ft 3 



(58) 



Here (• ■ •) denotes the average over noise realization. For 
simplicity, we assume that the fluctuating bias current is 
described as a Gaussian noise with the correlation func- 
tion of Eq. <|2fi|l and the spectral density of Eq. l(T7|) . We 
represent the average (• • •) by writing it as a functional 
integration over the noise. The transition probability V 
between different noise realizations may be described by 
the Fokker-Planck equation for the Ornstein-Uhlenbeck 
process^ 



dV 
~dt 



l_d_ 

Tndf 



(/>) 



dp 



(59) 



For this process, the transition probability V(f,t;f',t') 
for the noise from the value / at time t to the value /' 
after a time St = t — t' is given by 



V{f,t-J',t') 



2nnt 1 



-28tJT r . 



exp 



1 If - fe- dt/T "} 



-1/2 



St/r n ]2 



2n 2 



1 



~28t/r 



(60) 



We use this transition probability to express the proba- 
bility of specific noise realization as 

Vo(fo)V(fo,0;fut 1 )V(f 1 ,t 1 ;f 2 ,t 2 )---V(f n - U t n _ 1 ;f u t) 

(61) 
where f t = f(t t ) and V ? (f) = pTrnp" 1 / 2 exp(-/ 2 /2n 2 _) 
is the stationary Gaussian probability distribution of /. 
We note that 6U = (U — ij_i)/n. In the limit of SU — > 
(i.e., n — > oo), the average over the noise realization may 
be expressed as 



/ 1 V -L 



V 27m o 

(/ 2 +/ t 2 > 

e ~^~ 



<— I df Q df t D[f(t')} 



(62) 



where D[f] = nr=iW*/[47m 2 sinh^/^)] 172 } denotes 
the measure. The functional integral of Eqs. (|58|l and 
(|62|l is similar to that for a driven harmonic oscillator. 33 
Using this similarity, we may carry out the average over 
the realization of fluctuating bias current straightfor- 
wardly and obtain 



I(£)=I (f)expi -h, 

where i — t/r n , 
I (t) = e 



u>t 



coth 



(¥) 



(63) 
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(64) 
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Here the two dimensionless parameters, n — n /A and 
f n = t„A, characterize the amplitude and correlation 
time for the fluctuating bias current, respectively. For 
i *C 1, the fluctuation appears static. Hence the average 
(•••), which is over the static distribution of noises, yields 



7(t) = 



1 + i 



1 + t + i2n 2 T n t 



exp 



2tt 2 



(65) 



For i > 1, on the other hand, the fluctuation appears 
to be (5-function correlated (i.e., white noise), and hence, 
we obtain that I(t) oc exp(— t/T? mse ). This exponential 
decay indicates that the dephasing time TY mse may be 
expressed as 
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(66) 



We note that TV; olse is independent of T. Consequently, 
this will eventually cut off the divergent T? due to small 
coupling between the qubit system and environment at 
ultra-low T. Also, the effects due to the bias current is 
reduced, as expected, when asymmetry in the double-well 
potential vanishes (i.e., e = 0). 



V. DISCUSSION 

We now estimate numerically the coherence time for 
the JVQ and show that ultra-long coherence time can be 
obtained by using the experimental value for the param- 
eters. Since a long Nb-AlO^-Nb junction may be used 
to fabricate the qubit, we use the following experimental 
values in estimating T 2 ' 35 ' 36 Xl ~90 nm, A,/ ^25 /im, 
J c ~2xl0 6 A/m 2 and lo p ~90 GHz. For defmiteness, we 
chose a narrow width (i.e. L y ~0.2 /im) for the junction 
so that the quantum effect is enhanced. Also, we set that 
e/A = 0.01 since the variation in the pinning strength 
(i.e., e) for the two defect sites can be made small. Here 
we use the tunneling rate A which is obtained from quan- 
tum calculation of the ground state splitting, as shown 
in Fig. 4. 

As the coherence time T 2 depends strongly on A, we 
discuss, first, the dependence of A on the separation dis- 
tance £ between two defect sites and the pinning strength 
e. In Fig. 7, the tunneling rate for the ground state of the 
symmetric double- well potential (i.e., e = 0) is plotted as 
a function oi £ for e =0.21 (dashed line), 0.27 (solid line), 
and 0.33 (dot-dashed line). The curves show that MQT 
in LJJ depends strongly on both £ and e, as indicated 
by earlier studies^ The decrease in A with increasing i 
and/or e reflects that the tunneling rate decreases with 



FIG. 7: The dependance of the tunneling rate on the de- 
fect separation distance £ and on the pinning strength e is 
illustrated. 



increasing barrier potential V . We use this numerical 
result, below, in estimating T^. 

Ohmic environment: The relaxation time T\ and the 
dephasing time T? due to the interaction between the 
qubit system and environment are estimated by using 
Eqs. (|53|l and H56|) . respectively. These two characteristic 
times depend strongly on the quasiparticle dissipation ef- 
fect (i.e., j3). For defmiteness, we chose e = 0.27 and T — 
25 mK. This ultra-low T is chosen because the experi- 
ments show^ that the localized fluxon behaves as a quan- 
tum particle. For (£ — £ )/£ — 1-0, the estimated values 
are T x « (0.005//?) ns and Tf w (28//3) ns, indicating 

that both T\ and T^ become divergently long since /3 
is strongly reduced at this temperature. The estimated 
valued for j3 is roughly 0.03 at T ~ 4 K, and it is found 
to decrease exponentially with T^ below the supercon- 
ducting transition temperature. Phenomenologically, 38 
the dissipation effect represents the losses due to the tun- 
nel barriers. These losses are related to the quasiparticle 



resistance R qp {T) as 







UJpC R qp (T) 



(67) 



where C is the capacitance associated with the tunnel 
barrier. The T dependence of the quasiparticle resistance 
below the superconducting gap energy A sc (T) is given by 



R qp (T) = R T e 



A ac (T)/T 



(68) 



where Rt is the normal state tunneling resistance. The 
exponential T dependence for (3, as indicated in Eq. i|67|) 
has been observed to low temperatures^ (i.e., T <C T c ). 
We note that the dissipation coefficient (3 S , which rep- 
resents the contribution from the quasiparticle current 
along the junction layer, also decreases exponentially 
with T , 39 but this contribution is not included in this 
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work. This exponential T dependence for (3 suggests that 
both T\ and T? become divergently long at T — 25 mK 
since T x -v and T^ -> as T -> 0. Both T x and T^ 
are cut off by the decoherence due to weak bias current 
fluctuations. This suggests that the measured coherence 
time T 2 at T = 25 mK may be estimated as T 2 w T" olse 
since 1/T 2 ~ i/T" mse . 

Fluctuating weak bias current: The dcphasing time T^ 
is limited by the contribution due to the fluctuating bias 
current (i.e., current noise). As indicated in Eq. (|66H . 
the dephasing time due to the bias current fluctuation 
depends on the spectral density S no i se (0) = 2An 2 f n of 
Eq. I|27[) . The value of the spectral density S no i se (0) 
may be estimated by using the line- width, Sujppo, data 
for the Nb-AlO^-Nb flux flow oscillator (FFO). In esti- 
mating S no i se (0) we may use the relation between the 
magnetic field and bias current fluctuations. Recent 
measurements^ - of the line-width in FFO indicate that 
fluctuating bias current Sf(t) in the control line generates 
magnetic field fluctuation 5B(t) in the LJJ. In reverse, 
magnetic field fluctuations due to both external and in- 
ternal sources produce bias current fluctuations 4^ This 
suggests that, when magnetic field in the LJJ fluctuates, 
dephasing due to the bias current fluctuation may arise. 
Since the fluctuations are wide-band noises and are small, 
the relation between the bias current and magnetic field 
fluctuations may be expressed as 



SB{t) = K 5 f(t) 



(69) 



where K is the parameter of the order unity, describing 
conversion between the bias current and magnetic field 
fluctuations. We note that K depends on the geometry 
of the LJJ. The relation between the line- width Sujffo 
and the spectral density S no i se (0) for the current noise is 
given^i by 

2ir 
Sw FF „ = — (Rb + KR H ) 2 S noise {0) , (70) 
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where $ is the flux quantum in the superconducting 
state, Rb is the differential resistance associated with the 
bias current, and Rh is the differential resistance associ- 
ated with the magnetic field. The current noise spectral 
density measured from the FFO is S%£? e (0) ~ 2.8 x 1(T 22 
C 2 / s. Accounting for the geometry of the LJJ used in 
the experiment, the spectral density S^^ e (0) is given 
by d° e (0) = (2.5 x 10- 19 C 2 /s)S noise (0), Yielding 
S n oise{0) — 0.0011. Here we used the following param- 
eters that are obtained from the experimental data for 
the FF04& R B ~ 0.03 ft, R H ~ 0.005 ft, and K ~ 1. 
Also we chose Sujffo — 500 kHz, for definiteness, since 
the data indicate that the FFO line- width at the plateau 
of Fiske steps does not decrease below few hundred kHz 
at low values of Rb ■ 

In Fig. 8, we plot the dephasing time TV; olse versus 
the defect separation distance £ for e = 0.21 (dashed 
line), 0.27 (solid line), and 0.33 (dot-dashed line). Here 
rpnoise j g com p U ted from Eq. (|rjfj|> by using the experi- 
mental value of SVioise(O) — 0.0011 and by assuming, for 



definiteness, that r n = 1000 (Fig. 8a) and r„ = 500 
(Fig. 8b), which correspond to 10 ns and 5 ns, respec- 
tively. The curves indicate that TV; olse increases with r n . 
Also, T2° lse decreases and increases with £ and e, respec- 
tively, but it depends strongly on £ and weakly on e. The 
decrease in T^ mse with £ is due to the decrease in the 
tunneling rate, as shown in Fig. 7. We note that for 
(£ — £ )/£o = 1.3 the condition of f n ^> 1 is becoming dif- 
ficult to satisfy due to small tunneling rate. For r„ = 500 
and {£ — £ )/£ — 1-3, the computed value for f„ is 5.8, 
4.2, and 3.1 for e = 0.33, 0.27, and 0.21, respectively. For 
a smaller defect separation distance, say (l—t o )/t = 1-0, 
the condition of f n ^S> 1 is more easily satisfied and the 
computed value for T2° lse is in the microsecond range 
for both r n = 1000 and 500. For e = 0.27, T™ lse is 
roughly 55 [is and 30 /is for r„ = 1000 and 500, respec- 
tively. We compare this result with T^ for other super- 
conducting qubits. The measured values of T^ are 20 ns 
for the flux qubit^ at 25 mK, 10 ns for the phase qubit 3 
at 25 mK, and 500 ns for the quantronium 4 at 15 mK. 
These values indicate that T2 mse for the JVQ is orders 
of magnitude larger than the observed dephasing time 
in other qubits. This difference represents the fact that 
the fluxon's coupling to noisy environment is substan- 
tially weaker than other superconducting qubits. This 
indicates that longer coherence time may be obtained as 
fluctuating magnetic field in the LJJ is further reduced. 
Moreover, phcnomcnological comparison of these qubits 
indicates^ that the ultra-long dephasing time may also 
be attributed^ to the fact that the JVQ has much larger 
junction area than other qubits. 



VI. SUMMARY AND CONCLUSION 

In summary, we investigated the coherence time for the 
JVQ which may be fabricated by using a long Nb-AlO^- 
Nb juntion. Since the critical current fluctuation does 
not contribute to dephasing of the JVQ system, we esti- 
mate the coherence time by accounting for two sources of 
decoherence: i) quasiparticle dissipation and ii) current 
noise in the junction. We note that, within the low- 
est order approximation, the low frequency noise due to 
critical current fluctuation does not couple to the JVQ, 
and consequently it does not contribute to dechoerence. 
However, the low frequency noise due to bias current fluc- 
tuation is an important decoherence source. We showed 
that T\ and T? due to the quasiparticle dissipation (i.e., 
Ohmic environment) diverge at ultra-low temperatures 
(i.e., ~ 25 mK) since the dissipation effect (i.e., /?) be- 
comes exponentially small for T below the superconduct- 
ing transition temperature. In this case, the coherence 
time T 2 is determined by the bias current noise in the 
junction, as in many superconducting qubits. We esti- 
mated X™ se by accounting for the fact that the current 
noise may arise from the magnetic field fluctuations in 
the junction. This bias current fluctuation is described 
realistically by using the Gaussian colored noise with a 
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FIG. 8: The dephasing time T^ mae versus the defect sep- 
aration distance £ is plotted for e=0.21 (dashed line), 0.27 
(solid line), and 0.33 (dot-dashed line) to illustrate the de- 
pendence of T% mse on both t and e. The computed T%° ise 
for a) t„ — 1000 and b) 500 illustrate the dependence on the 
correlation time. 



long correlation time. Our estimated value of TV; mse for 
the JVQ, which is obtained by using the experimental 
data from the Nb-AlO^-Nb FFO, is in the microsecond 
range because the spectral density for fluctuating mag- 
netic field is very low. The value for T<p is few orders of 
magnitude larger than that measured for the quantron- 
ium, suggesting that the JVQ may also be a good can- 
didate for quantum computer. This surprisingly long co- 
herence time for the JVQ is due to the fact that the 
fluxon, which behaves as a topologically stable quantum 
particle at ultra-low temperature, couples very weakly to 
both internal and external noise sources. The current es- 
timate for Ti may be extended further if the magnetic 
field fluctuations, which are considered as the dominant 
decoherence source in LJJ, can be further reduced. 



This work suggests possibility that a superonducting 
qubit with an ultra-long coherence time may be realized 
by exploiting quantum property of fluxon pinned in a 
double-well potential in LJJ. This work also provides in- 
sight into design and fabrication of the JVQ. The cur- 
rent approach may be easily extended^ to realize mul- 
tiple noninteracting qubits in quasi-one dimensional LJJ 
stacks. Hence, it would be interesting to verify macro- 
scopic quantum coherence behavior by either spectro- 
scopic measurement of level splitting or by observation 
of Rabi oscillations in the JVQ. 
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